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Abstract 

Recent experiments in a mixture of two hyperfine states of trapped Bose gases show behavior 
analogous to a spin-1/2 system, including transverse spin waves and other familiar Leggett-Rice- 
type effects. We have derived the kinetic equations applicable to these systems, including the spin 
dependence of interparticle interactions in the collision integral, and have solved for spin-wave 
frequencies and longitudinal and transverse diffusion constants in the Boltzmann limit. We find 
that, while the transverse and longitudinal collision times for trapped Fermi gases are identical, the 
Bose gas shows diffusion anisotropy. Moreover, the lack of spin isotropy in the interactions leads 
to the non-conservation of transverse spin, which in turn has novel effects on the hydrodynamic 
modes. PACS numbers: 03.75.Mn,05.30Jp,05.60.Gg,51.10.+y,67.20.+k. 



1 



In recent JILA experiments, a mixture of two hyperfine states was found to segregate 
by species. The theoretical explanationj^-j^ for this behavior is based on the two states 
playing the role of a pseudo-spin- 1/2 system, having transverse spin waves. The theory of 
these new effects is based on old ideas of the transport properties of polarized homogeneous 
quantum gases of real spins, such as ^He gas and solutions of ^He in liquid ^He-lgj]. jl^ 
transcribed to the trapped gas pseudo-spin case. 

Besides spin waves, the theory for homogeneous polarized fermions or bosons led to 
the prediction of anisotropic spin diffusion in the degenerate state.P],[l^,|lll When a spin 
nonuniformity is longitudinal, that is, with a variation in the magnitude of the magneti- 
zation, the spin diffusion coefficient is D^\. On the other hand, in a spin-echo experiment, 
the magnitude of the magnetization is uniform but it varies spatially in direction. The 
corresponding diffusion coefficient, D± is less than D\\ when the system is polarized and 
degenerate. Experimentally this feature has been seen, but was not always in reasonable ac- 
cord with theory, Moreover, Fominjl^ has suggested the effect should not exist. However, 
a recent experiment!^ has overcome several possible experimental objections and finds good 
agreement with theory. Moreover Mineev has very recently presented theoretical analysis 
that questions the validity of Fomin's approach. 13| 

Thus it seems useful to see whether a similar difference between longitudinal and trans- 
verse diffusion in trapped gases might provide an alternative testing ground for this question. 
However, what we show here is that the physical possibility of having differing interac- 
tion parameters between up-up, down-down, and up-down states (interaction anisotropy) 
provides a new physical basis for anisotropic spin diffusion for bosons even in the Boltz- 
mann limit. For longitudinal diffusion in the Boltzmann limit only up-down scattering 
contributes. However, in the transverse case, two spins at differing angles approach one 
another, and the scattering can be analyzed as being a superposition of, say, up-up and up- 
down scattering. In the fermion s-wave case, the up-up part gives no contribution, and, in 
the Boltzmann limit, the diffusion coefficients are identical. In that case one must go to the 
degenerate limit to see the anisotropy, which then is expected to arise because the density 
of scattering states differs in longitudinal and transverse cases. On the other hand, for 
bosons, for which both the up-up and down-down scattering rates do contribute, we find an 
anisotropy even in the Boltzmann limit, but only if the various scattering lengths differ. We 
have here the striking effect that, although both gases obey Boltzmann statistics, there is a 
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macroscopic difference between fermion and boson behavior. 

The presence of interaction anisotropy provides another unusual effect, namely that trans- 
verse spin is not conserved. j3| This leads to a decay of the transverse spin (a T2 process) that 
seriously affects the hydrodynamic modes of the system. Below we first use the moments 
method to compute the spectra of the lowest-lying longitudinal and transverse modes. How- 
ever, with that method we obtain a transverse decay rate 7^ that diverges as r approaches 
zero, in contrast to the usual diffusive behavior where 7^ oc r. In this case it is necessary 
to solve the local hydrodynamic equations to find the correct behavior, in which the hydro- 
dynamic solutions are localized at the low-density regions at the edges of the cloud where 
the collision time is longer. The result is a much smaller decay rate than that obtained with 
the moments method. 

In our previous work, Ref. we derived an analog of the Landau-Silin equation for a 2 x 2 
density operator rip (here acting in the pseudo-spin space), with effective mean- field single 
particle energy matrix ip. We can write the density and single-particle energy in a Pauli 
representation as rip = ^ (^fpl + lUp ■ and ip = {cpl + hp ■ ct^ where 5" is a Pauli matrix, 
|(/p ± rripz) give the diagonal components of the density ripi = ripu, while nip represents the 
polarization, which in equilibrium is along the axis z. We find the following approximate 
equation for lUp-. 

Pi duip dU dm.p 



dvcvp 2 ^ 
^-^hxmp + E 



m dvi dri dpi 



Tr {aip} (1) 



with mpzi^r) = ripi — np2 and Upuir) = n*p2i = |mp_(r) = |(mp^. — irripy). The 2x2 colhsion 
integral is Ip. The effective mean magnetic field 

h = ^z + r/yM (2) 
where hQ.o = Vi-V2 + [(tn - t^) rii - (^22 - ^12) ^^2] 

X (l-|-?7). In these rj is 1 (—1) for bosons (fermions); M(r) = / dp/h^ mp(r); ni(r) = 
/ dp/h^ np{r)] Vi is the external field for species i; U = |(Vi + V2); and Mz = ni — n2. The t's 
can be evaluated in terms of the measured scattering lengths a^/? by using ta/3 = ^Trhaa/s/m. 

The equilibrium solution in the Boltzmann limit is m^^ = A4(PhujY exp[— /5(j9^/2m-|-f/)] 
where N is the total number of particles, Ni is the number of species i, Ai = Ni — N2 is the 
total magnetization, and cD = (tOx^ytOzY^^ . 
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We have derived the colhsion inteErral for the Boltzmann case when the various interaction 
pa..a.ente. diffe. „o. ag.ee. with the sa.e quantity derived in Refs.QandB 
and reduces properly to previous results if all the t's are taken equal. [3| , [isl We find 



= ^ / dpidp2dp35{pi + P2 - P3 - P4)(5(epi + Cpai - ep:n - ep4y 



P2j<J2Cr2 



[(^Pl)o-o-'('^p2)o-20-2 +^('^Pl, 



(JCT2 l^p2 /0-20-', 



~\~'^taa2^u'a2 [(^P3 )o-cr' ('^p4 )(T2a-2 ~l~ ''^('^ 



P3y'o-(T2 V^P4y'o-2 



(3) 



We will linearize the kinetic equation for nip around the global equilibrium value m^^z 
and use a moment approach to compute the spin wave and diffusive damping just as done 
previously. As in Ref. |5| we assume that the effective longitudinal field Q,q can be 

adjusted experimentally to zero. The linearized longitudinal and transverse equations are 



dSm. 



pe 



dt 



+ E 



Pi d5m. 



pe 



dU d5m 



pe 



m dvi dvi dpi 



Y,(^{(^\Lp\cr) 



(4) 



and 



d6m. 



'P+ 



dt 



+ E 



Pi d6mp^ dU d6m 



p+ 



m dvi dvi dpi 



2(2|4|1) 



(5) 



where Lp is the linearized form of Ip 



In the following, for brevity, we compute only results for the monopole and dipole modes 
although experiments have detected the quadrupole modes. Similar arguments hold for the 
quadrupole case, which we will present in a longer publication. 

Longitudinal case: We use a variational function of the form 

Srripz = (do + ciiZ+a2Pz)m'p^ (6) 

and take the 1, z, and pz moments of the kinetic equation in both the longitudinal and 
transverse cases. The results for the longitudinal case, if we assume a time dependence of 
exp{iLjt) for ai and 02, are 



daQ I dt 
iooai — 002(^2 
i(jja2 + oOztti 
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-711^2 



(7) 
(8) 
(9) 



with 7|| = 470/3 where 70 = 7if3m^u^ti2N/h^ comes from integrating the colhsion integraL 
Eq. (|7j) indicates that the monopole mode does not decay in the longitudinal case, which 
is consistent with the conservation of magnetization. The second line is the magnetization 
equation of continuity. The relaxation rate 7|| agrees with that derived in Ref. Q. The dipole 
spectrum is plotted in Fig. ^ as a function of = l/7j|, the spatially averaged collision time. 
In the small t^^ limit, one finds 

W||=iwfr||, (10) 

which has the form of the lowest-order solution of a diffusion equation in a harmonic poten- 
tial. 
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FIG. 1: Real (dashed) and imaginary (solid) components of longitudinal dipole spin wave modes 
versus average relaxation time ry. Note the linear dependence of Im(u;) for small ry characteristic 
of diffusive behavior. The dash-dotted lines represent the results of a numerical calculation to be 
discussed below. 

Transverse case: We again use the form of Eq. Taking 1, z, and Pz moments of 

Eq. (0) yields the results 

dao/dt = — 7Tao (11) 

1 , , 

iLjai-uJza2 = --'Jrai (12) 
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— Co'Af)«2 + t^2«l = — 7±a2 (13) 

where 7t = 7o(l + v) {^^i^r^f U with U = NJN, and 



7± = 7|| 



7R - 3S 



(14) 



^12 

with R = {l+r]) J2cTtlafa + {l-v)tl2 and S = 2ti2[{l+ri) J2ataafa-Vti2], and the mean-field 
frequency is 

Its j '''' 

where A is the thermal wavelength. 
Comments: 

1) If the interactions parameters tij are all equal, we have 7t = 0, i? = S* = so that 
7j_ = 7||. Eqs. (fTT| - (fT^ then reduce to those of Ref. and the longitudinal and transverse 
relaxation rates are the same, which agrees with the standard result for a homogeneous real 
spin system in the Boltzmann limit. 

2) For fermions, we have t] = —1, so that, even if the t's are not equal, 7t = and 
7± = 7||- 

3) For bosons with unequal t's, the spatial averaged transverse relaxation rate is not 
generally the same as the longitudinal. Moreover, we have a T2-type relaxation rate for 
ao and in the equation of continuity ()12j) . The interaction anisotropy behaves something 
like a dipole-dipo 
previously in Ref. 

If, for now, we take jt = 0, then the lowest mode in the hydrodynamic limit takes the 
form 

where = 1/7_l and the so-called "spin-rotation parameter" fi = uj mT \ . The form of 



e interaction allowing relaxation of the transverse spin, an effect noted 

3. 



Eq. (|T6|l is the hydrodynamic frequency as modified by spin rotation. P , (10] , [l^ • The first 
term is the effective diffusion frequency while the second is the dipole-mode pseudo-spin-wave 
frequency. 

The effect of non-zero 7r is to allow a T2 relaxation of the transverse spins. The results 
are shown in Fig. 2, where we have set I/t^ = 7||, 7r = 0.02/r_|_, and uom = '^z- In the small 
limit, one no longer has the hydrodynamic decay rate approaching zero, but instead it 
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FIG. 2: Real (dashed) and imaginary (solid) components of transverse dipole spin wave modes 
versus average relaxation time t±_ when the transverse spin is not conserved. The dotted line 
shows the lower imaginary mode when the transverse decay rate 7t = 0. The mean-field frequency 
ujm is taken as ojz- The linear behavior of Im(a;) for small t± characteristic of diffusive behavior is 
destroyed and replaced by a divergence within the moments method used here. 

diverges at the origin because Im(a;) (o^^tj, + 7t) , and 7t ~ 1/tj_. However, although 
suitable for finite lot, the moments method is inadequate in the hydrodynamic limit. It 
fails because the simple forms assumed for spatial dependence cannot adjust to the spatially 
dependent relaxation rates. One must solve local equations numerically for the spatial 
behavior. 

To obtain the hydrodynamic equations we expand the momentum distribution in terms 
of Hermite polynomials 



fc=0 

Substituting this into the kinetic equations, integrating over the momentum, and keeping 
terms lowest order in t±_ gives, in the transverse case. 




(17) 



dt5M+ + d,J+ = -jt{z)SM+ 
kT 

dtJ+^ dJM+ + ujlz5M+ + iujM{z)J- 



(18) 
(19) 



where SM+{z,t) = co{z,t) is the nonequihbrium magnetization density, J{z,t) = 
J dp/h^{p/m)6mpz = ci{z,t) is the spin current, and7_L(2;) = 'y±{0) exp^—Pmuj^z^ /2). Anal- 
ogous equations hold in the longitudinal case. On the RHS of Eq. (jl9p the k = 1 momentum 
distribution has been treated as an eigenfunction of the linearized collision integral. This is 
justified by a numerical calculation of the matrix elements of the collision integral, which 
gives 

L^[H,{p)] = -7^(0)(1.000/7i(p) + 0.123/73(p) 

-0m09AH^{p) + ...) ^ -7x(0)ifi(p) (20) 

The eigenvalues of the hydrodynamic equations have been calculated numerically for the 
dipole mode with boundary conditions 6M{0) = 0, J(0) = 1, and J(oo) = 0, and the 
monopole mode with boundary conditions 6M{0) = 1, J(0) = 0, and J(oo) = 0. For the 
longitudinal and isotropic transverse cases this leads to only small corrections to the r — > 
part of the spectra obtained by the moments method as shown in Fig. [TJ 
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FIG. 3: Imaginary part of the spin-wave spectrum vs. uJzT± for the monopole and dipole modes 
with lom = ^^z, 7t = 0.027^, and 7± ~ 7|| for both the moments method (thick) and hydrodynamic 
calculations (dashed). 

However, for 7^ > the hydrodynamic spectrum differs qualitatively from that of the 
moments method calculation. As — the hydrodynamic dipole and monopole modes 
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FIG. 4: Profiles |(5M(z)| of the the dipole modes vs. y/ (3mu-^z, normahzed with a Hermite weighting 
function. Tallest peak is the hydrodynamic mode for 77^ = and t± = 0.01. Middle peak is the 
moments method ansatz. Smallest peak is the hydrodynamic mode for jt = 0.027_|_ and w^r = 0.01. 

do not decay at a rate ~ but instead decay at a slower rate ~ y'log {1/u!zT±) (See 

Fig. El) In fact, at small enough uJzT± the T2 decay of the magnetization at the center of 
the trap causes the monopole and dipole modes to coalesce into spin-waves localized on the 
lower density regions on the left and right sides of the trap. (See Fig. EJ) 

In experiments on Rb, the interaction anisotropy is very small. To test the novel effects 
predicted here one might use Na, which has a difference in interaction paraments; Numer- 
ically we estimate that for ^^Na 7^ can differ from 7^ by as much as 14% with 7t/7± ~ 0.04. 
Interaction differences might also be induced by using Feshbach resonance methods. 

We thank Dr. Jean-Noel Fuchs and Prof. David Hall for useful discussions. 
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